Charged Pion Production in 2 to 8 AGeV Central Au+Au Collisions 
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Momentum spectra of charged pions over nearly full rapidity coverage from target to beam rapidity 
have been extracted from 0-5% most central Au+Au collisions in the beam energy range from 2 
to 8 AGeV by the E895 Experiment. Using a thermal parameterization to fit the transverse mass 
spectra, rapidity density distributions are extracted. The observed spectra are compared with 
predictions from the RQMD v2.3 cascade model and also to a thermal model including longitudinal 
flow. The total 47r yields of the charged pions are used to infer an initial state entropy produced in 
the collisions. 

PACS numbers: 25.75.-q, 25.60.Gc 



I. INTRODUCTION 

One of the primary goals of the study of heavy ion col- 
lisions at relativistic energies is the improvement of our 
understanding of the bulk properties of nuclear/hadronic 
matter at high temperatures and densities. As a first step 
in understanding these properties, one should carefully 
characterize the particle species that make up the bulk 
of the matter. For the energy regime of the Bcvalac/SIS 
(0.2 to 1.2 AGeV), collisions between heavy nuclei cause 
a compression of the nuclear matter, resulting in a disas- 
sembly into the constituent neutrons and protons, which 
are emitted either individually or bound within various 
light composite fragments (d, t, 3 He, 4 He)0,0- 

At the top energy of AGS Au+Au collisions of 10.8 
AGeV, the most copiously produced charged particles are 
the lightest of the mesons, the pions0, 3- The matter 
has evolved from a heated and compressed gas of nucleons 
into a hot dense gas of hadrons, predominantly pi mesons. 
Thus the energy region from 2 to 8 AGeV represents a 
transition regime. Studies of nuclear stopping suggest 
the maximum density achieved increases from three times 
normal nuclear density at 1 AGeVyjJfjl, |jfl to eight times 
normal nuclear density at 10 AGeV[8j. Measurements of 
the proton directed and elliptic flow indicate that hydro- 



dynamic flow has saturated across this energy range and 
cannot account for the increased energy available |& llPf. 
This additional available energy goes primarily into pion 
production. By studying the pion production across this 
regime we are able to observe nuclear matter in transi- 
tion. 

Early on in the study of heavy ion collisions, an en- 
hanced yield of pions was seen as a possible signature 
of a transition to a deconfined state of matter. How- 
ever, in the early studies of pion yields in the 1 AGeV 
energy range at the Bevalac^J Il2l Il3j . it was ob- 
served that the measured pion production cross sections 
were smaller than predicted at the time. This obser- 
vation led to the conjecture, which was later experi- 
mentally demonstrated, that the excess kinetic energy 
was converted into hydrodynamical flow effects. The 
strong radial flow observed at this energy implies a sig- 
nificant expansion and cooling, which limits the freeze- 
out pion multiplicities[T ^. Il5 | . More recently, more de- 
tailed mesurements of the pion yields from 1 AGeV 
Au+Au collisions have become available from the SIS 
experiments [Hi ITtI ll8L ITgj . These results demonstrate a 
roughly two to one ratio of 7r _ over 7r + and a strong non- 
thermal low-p t enhancement. Both features are strong 
indicators that pions at this energy are produced almost 
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exclusively through the Delta resonance [20j. In full en- 
ergy AGS collisions (Au+Au at 10.8 AgeV), pion pro- 
duction has been studied at midrapidityQ and at target 
rapidity |2l|. Although there is still an asymmetry be- 
tween 7r~ and 7r + production and there is also still evi- 
dence of a low-pt enhancement in the momentum spectra, 
these effects are much less significant than at SIS. The 
broad rapidity coverage for both pions and protons at 
these energies has also been used to study the develop- 
ment of longitudinal flow[22j. 

This paper will detail the development of pion produc- 
tion across the beam energy range from 2 to 8 AGeV. The 
role of the Delta resonance production mechanism will be 
explored through observations of the overall pion ratios 
and through RQMD simulations. The rapidity density 
distributions will be used with previously published pro- 
ton rapidity distributions |2^| to explore the effects of col- 
lective longitudinal flow. The overall pion multiplicities 
will be used to establish a low energy baseline for a re- 
cent pion multiplicity based QGP search at the CERN 
SPSS- 



II. DATA COLLECTION AND ANALYSIS 

The data were taken at the Brookhaven National Lab- 
oratory Alternating Gradient Synchrotron (AGS) by the 
E895 Experiment using the EOS TPC [23 during a series 
of runs in 1996. This article presents charged pion trans- 
verse mass and rapidity density spectra from Au+Au 
collisions at nominal beam energies of 2, 4, 6, and 8 
GeV/nucleon (AGeV). (After correcting for energy loss 
before the target, the actual beam energies of collisions 
at 2 and 4 AGeV were found to be 1.85 and 3.91 AGeV, 
respectively. No corrections were necessary for 6 and 
8 AGeV.) The EOS TPC provides nearly 4tt solid an- 
gle coverage, which makes global characterization of the 
collision events possible. Charged particle momenta are 
reconstructed from the helical trajectories of tracks re- 
constructed from the ionization trails left by particles 
passing through the TPC, which was situated inside the 
Multi-Particle Spectrometer (MPS) magnet. Data pre- 
sented for 2 AGeV collisions were taken in a 0.75 Tesla 
field, while 4, 6, and 8 AGeV collisions were taken in a 1 
Tesla field. 

A primary track multiplicity for each event is obtained 
by rejecting those tracks which do not pass within 2.5 
cm of the reconstructed event vertex in the target. The 
multiplicity distribution from a minimum trigger-biased 
sample of events is used to discriminate event centrality 
classes by assuming a monotonic relationship between the 
impact parameter and multiplicity |26|. The data pre- 
sented in this paper are selected on the top 5% most cen- 
tral events; an ensemble of approximately 20,000 events 
at each energy for this centrality selection were used to 
obtain the spectra. 

Identification of particle species is determined via mul- 
tiple sampling (up to 128 samples) of the ionization in 



P10 (10% methane, 90% argon) drift gas. The average 
ionization energy loss, (dE/dx), is computed for each 
track from the available samples via a truncated mean 
to reduce the influence of the large energy tail of the dis- 
tribution. The 20% highest dE/dx samples are discarded 
from the calculation of the average. 

Fig. shows a typical particle identification map with 
(dE/dx) plotted versus reconstructed track rigidity (= 
p/Z) at 6 AGeV. The various particle species are iden- 
tified by their separation into bands. The pion spectra 
were obtained by fitting projections of the (dE/dx) in 
narrow bins of mt — mo and rapidity, using an assumption 
of the pion mass and charge to calculate the (m t — mo,y) 
coordinates for a given measured particle (p x ,p y ,p z ). 

The single particle (dE/dx) projections are often de- 
scribed by a Gaussian distribution centered on the mean 
value predicted by a Bethe-Bloch formulation. However, 
this assumes that the calculation of (dE/dx) for each 
track was obtained in an identical fashion for all tracks 
in the distribution. In fact, the truncated mean method 
used in this analysis introduces a skewing toward larger 
(dE/dx) which comes from combining tracks of differ- 
ent number of samples, ^hits- One way to avoid this 
skewing is to divide the data into bins of N^i ts , which 
reduces the effect. However the reduction in statistics 
for each bin leads to increased uncertainty in the deter- 
mined yields. Therefore, for this analysis, a model of 
the N/jits-integrated single particle (dE/dx) distribution 
shapes was used with the predicted mean values from a 
Bethe-Bloch parameterization of the (dE/dx) as a func- 
tion of /?7 to extract the total yields of pions from each 
projection. The model is represented by a correlated sum 
of Gaussian distributions: a main Gaussian for the bulk 
of the distribution plus a smaller, offset "shoulder" Gaus- 
sian for the high-(dE/dx) tail. The parameters of the 
model were studied as a function of beam energy, parti- 
cle type and (m t — mo,y) bin in regions where the sin- 
gle particle distributions can be separately characterized. 
Tight constraints on the model parameters were applied 
in order to extrapolate the model into regions where mul- 
tiple particle distributions partially overlap. Thus the 
relative yields of different particles were deconvoluted in 
each (mt — mo,y) bin: 0.1 unit rapidity slices over the 
full rapidity range from target to beam rapidity, with 
the midrapidity bin covering the range -0.05 < y cm < 
0.05, and in 25 MeV/c 2 bins in m t -mo in the range 
< m t -mo < 1.0 GeV/c 2 . Particles of the wrong mass 
and/or charge contaminating the pion sample in a given 
(mt — TOo,y) bin were discarded. 

Fig. |21 shows an example of a single (mt — mo,y) slice 
at midrapidity for 6 AGeV and 0.125 <m t — m$< 0.150 
GeV/c 2 . The total yields are obtained by integrating the 
area under the fitted distributions for each particle. Op- 
positely charged particles are projected separately. The 
inset shows the projection of negatively charged particles, 
which are well-characterized by this method of particle 
identification over the entire range of transverse mass and 
rapidity. However, as the rapidity and transverse mass 
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FIG. 1: (dE/dx) as a function of rigidity at 6 AGeV. A pa- 
rameterized Bethe-Bloch prediction is used to fix the mean 
energy loss for each particle species as a function of rigid- 
ity. The narrow slices illustrate the location of the (dE/dx) 
projections shown in Fig. 



increase, the positively charged pions suffer increasing 
contamination due to the overlap with the positive kaons 
and protons. The range of m t — too and rapidity over 
which the ir + yields are extracted is therefore more lim- 
ited. 

Electron contamination was estimated by studying 
their yields in the regions of phase space where they 
are relatively cleanly identified (for pi a b < 150 MeV/c 
and 300 < p; a 6 < 500 MeV/c) and interpolating between 
these limits. It was found that electrons contribute ap- 
proximately 10% to the observed yield of pions for lab 
momenta 150 < pi a b < 250 MeV/c. Since the electron 
yields fall significantly as a function of their own trans- 
verse momentum, the contamination predominantly ef- 
fects the lowest pion m t — too bins. The errors on the 
quoted pion yields account for possible systematic uncer- 
tainties in the determination of the electron contamina- 
tion. 

Observed kaon yields from the same beam energy 
range, measured by E866/E917 [Hill, folded with the 
EOS TPC detector response allows us to extend the reach 
of the 7T + fitting out to pi a b ~ 1.2 GeV/c. For lab mo- 
menta above 1.2 GeV/c, however, the ir + become hope- 
lessly entangled with the protons, which are approaching 
minimum ionizing (dE/dx). This cut-off manifests at 
different m t — too, depending on the rapidity and bom- 
barding energy. 

Detector response to the final-state collision products 



FIG. 2: Projections of the positive and negative (inset) parti- 
cle (dE/dx) at mid-rapidity and 0.125 < rat — mo < 0.150 
Gev/c 2 , in the rigidity ranges indicated on Fig. The 
(mt — rao,y) ranges were calculated assuming the pion mass. 
The function used to fit the data is the two-Gaussian model 
(described in the text) of the single particle (dE/dx) distri- 
butions, one for each particle type. The integrals under the 
pion peaks are the total yield for the fitted (m t — mo,y) bin 
and the contributions from all other particles are discarded. 



has been extensively studied using the GEANT 3.21 sim- 
ulation package. Small (maximum of 4 per event) sam- 
ples of pions with momentum distributions approximat- 
ing the real data are embedded into full data events. 
These particles are tagged and propagated through the 
data reconstruction chain to determine the effects of de- 
tector acceptance, tracking efficiency, and momentum 
resolution. 

Since the beam actually passes through the sensitive 
volume of the detector, there is no explicit low-pi mea- 
surement cut-off. However, forward-focusing causes in- 
creased tracking losses at low p t and forward rapidities 
due to higher track densities and track merging. Losses at 
backward rapidities and high m t — too are dominated by 
the geometric acceptance of the detector, while at more 
forward rapidities, the high m t — mo tracks are stiffer 
and therefore suffer from worsening momentum resolu- 
tion. Fig. shows the overall detection efficiency ob- 
tained from simulations as a function of m t — Too and 
rapidity for pions at each beam energy. The contours 
indicate steps of ~15%. The raw data are corrected for 
these effects to extract the total yields of positive and 
negative pions as a function of m t — too and rapidity at 
each beam energy. 

Due to the large acceptance of our device, we can test 
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III. RESULTS 
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FIG. 3: EOS TPC total "efficiency" (= 1/correction, includ- 
ing effects of acceptance, tracking efficiency and momentum 
resolution) for pions produced in 2, 4, 6 and 8 AGeV 0- 
5% central Au+Au collisions. The efficiencies are obtained 
from GEANT simulations of the detector response to Monte 
Carlo tracks embedded into real data events and propagated 
through the reconstruction chain. The contours, shown in 
step of ~15%, include the effects of geometric acceptance, 
tracking efficiency and momentum resolution. 



Figs. 0] and [5] show the fully corrected, invariant yields 
of charged pions per event from 0-5% central Au+Au col- 
lisions at 2,4,6 and 8 AGeV. Mid-rapidity is shown un- 
sealed as black circles, while each bin forward/backward 
of mid-rapidity is scaled down by a successive factor of 10. 
Forward rapidities are indicated as open circles and back- 
ward rapidities as open triangles. The reported error bars 
include both statistical and systematic uncertainties. 

The approximately exponential decay of the particle 
yields as a function of transverse mass has been observed 
in high energy particle and heavy ion experiments over 
a wide range of conditions. In order to extract the full 
4-7T yields of pions (including those from resonance feed- 
down) produced in the beam energy range studied here, 
a simple parameterization of the pion m ( — too spec- 
tra which reproduces the observed shapes of the spectra 
over the full range of m t — to , from to 1 GeV/c 2 is 
used. This parameterization is the sum of two indepen- 
dent Maxwcll-Boltzmann thermal functions, each term of 
which can be expressed in terms of the measured coordi- 
nates, (mt — TOo,y), integrated over azimuth as 
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d 2 N 



27TTOt dm t dy 



A{y)m t e^ i - mt - m ^/ T '"^ (1) 



where the amplitude A and inverse slope parameter, T e ff 
are parameters which can be extracted from fits to the 
data at each rapidity slice. Integration of either of these 
two distributions over m t produces its contribution to the 
total number of particles per unit of rapidity in the given 
rapidity slice [U- 



the systematic uncertainties on the yields by comparing 
forward and backward rapidity corrected spectra. Over- 
lapping acceptances at 2, 4, 6 and 8 GeV allow us to 
conclude that rapidity bins corresponding to yi a b < 0.5 
are most affected by the systematic uncertainties in our 
corrections. The underestimated corrections in this re- 
gion of phase space are not surprising given that nei- 
ther detector performance nor detector simulations were 
optimized for target rapidity in the lab frame. Here 
tracks cross the fewest pad rows and have the small- 
est radii of curvature. Mid-rapidity yields at all four 
beam energies were checked against published results 
from E866/E917[2i|33. The average level of agreement 
between our spectra and E866/E917 over all m t — too and 
beam energies is observed to be approximately 7%. The 
minimum systematic uncertainties at all rapidities for the 
E895 data presented here are estimated to be 5%, while 
for the most backward rapidity bins, the uncertainties at 
high to 4 — m can become as large as 50%. These errors 
are included on the spectra reported in this article. 



^(y) = 2nA(y)(mlT eff (y) + 2m T* ff (y) + 2T* fJ (y)) 

(2) 

Note that A(y) in Eq. (JTJ can be re-written in terms of 
dN/dy(y), using Eq. 0. 

The low-pi enhancement observed in the pion yields 
has been attributed to feed-down from late-stage reso- 
nance decays [32l H3> 03 j which tend to populate lower 
pion p t in the frame of the collision due to the nature 
of the decay kinematics. In particular, at the beam 
energies presented here, the delta resonances, (such as 
the A(1232)) are the predominant mechanism for pion 
production, due to the very large cross-section for pion- 
nucleon interactions. The observed asymmetry in the 
spectral shapes of positive and negative pions at low 
nit — wo has been described elsewhere by a final-state 
Coulomb interaction of the pions with the nuclear fire- 
ball mm mm. 

The two-slope model which has been applied to the 
pion spectra was chosen to provide the simplest phe- 
nomenological description of the data. Since this pa- 
rameterization does reasonably well at describing the ob- 
served shapes of the spectra over all m t — too, it is used 
to extract the in yields of pions in these collisions with 
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FIG. 4: Invariant yield per event as a function of m t — mo for n~ at 2, 4, 6, and 8 AGeV. Midrapidity is shown unsealed, while 
the 0.1 unit forward/backward rapidity slices are scaled down by successive factors of 10. The functions plotted with the data 
are two-slope Boltzmann parameterizations described in the text. 



minimal extrapolation. 



d 2 N _ dN/dy 1 m t e- ( - mt - m °'>/ Tl dN/dy 2 m t < 



2-Kmt dm t dy 2-kT x (rag + 2m T 1 + 2T{) 2-kT 2 {ml + '. 

The four parameter fit includes two independently fit 
inverse slope parameters (Ti(y), T2(y)), and two inde- 
pendently fit yield parameters (dN/dyi(y), dN/dy2(y)), 
which reasonably describe the low- (mt — mo) and high- 
(m t — mo) portions of the spectrum at each rapidity. 
These fits are shown as solid lines in Figs. 0] and [5] 
The data at more forward rapidities, where in partic- 
ular the 7r + spectra are limited to the range m t — mo < 
200 MeV/c, are not included in the two-slope parameter- 
ization fits. 

More detailed examples of these two-slope fits at mid- 
rapidity are shown in Fig- EI Figs. I7I1T3I show the four fit 
parameters as a function of rapidity from the n~ spectra 
and the resulting total dN/dy. The parameters plotted in 
the top row are from fits in which all four parameters are 
allowed to be free. All of them show an approximately 
Gaussian dependence on rapidity, which is demonstrated 
by the solid line in each panel. The bottom row of each 
figure demonstrates the result of a second fit of the spec- 
tra in which the high- (rat — mo) inverse slope parameter 



is constrained to the Gaussian (solid line) value. This 
procedure smooths out the covariance between the in- 
-(ntitvwlWfa yield parameters but has a negligible effect on 

2mo l 7 J 2 La r ; $ffif y ( ri S ntmost P ane1 )- 

If the source of pions were a static thermal source 
of temperature To and zero chemical potential, the ex- 
pected shape of T e ff in Figs. 1711 II would be To/ cosh(y). 
(Of course, there would then be needed only a single term 
in Eq. ©.) We have shown in Ref. that there is 
significant longitudinal flow - i.e., the source is not static. 
As the main purpose of the fits in Figs. 21and[S]is f° r in- 
tegration of the (m t — mo,y) spectra to obtain 47r yields, 
a Gaussian fit of the inverse slope parameters in Figs. 
171111 rather than l/cosh(y), was used. 

Fig. II II shows the ir + fit parameters at each beam en- 
ergy. The limited range in rat — mo of the tt + spectra is 
caused by the (dE/dx) overlap with the protons, which 
makes it impossible to get good spectral data at all rapid- 
ity slices. In order to make reasonable estimates of the 
4-7T positive pion yields, the high-(m t — mo) positive pion 
inverse slope parameters were assumed to be the same as 
those of the negative pionspHj. In the estimation of the 
7r + dN/dy, the high- (mt — mo) inverse slopes were there- 
fore fixed to the negative pion values. Experimentally, 
the mt — mo negative pion inverse slope parameters do 
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FIG. 5: Invariant yield per event as a function of m t — mo for ir + at 2, 4, 6, and 8 AGeV. Midrapidity is shown unsealed, while 
the 0.1 unit forward/backward rapidity slices are scaled down by successive factors of 10. The functions plotted with the data 
are two-slope Boltzmann parameterizations described in the text. 



reasonably reproduce the observed high- (mt — mo) posi- 
tive pion inverse slopes in the rapidity regions where they 
can be measured. Where the m t —m spectra are severely 
truncated, the positive pion yields were obtained from 
single-slope fits. The reported systematic uncertainties 
in these rapidity slices are correspondingly larger to ac- 
count for this missed yield. 

Fig. E| shows the beam energy dependence of the ra- 
pidity distributions of negative and positive pions ex- 
tracted from the fits. These distributions are well- 
described by Gaussians, which are used to obtain the 
total 47T yields by integrating the fitted distributions over 
all rapidity. 

(tt) = / — dy = V^wQ (4) 
J-oo dy 

where Q is the value of dN/dy at y-y cm = and w is 
the width. The fit parameters and their uncertainties are 
listed in Tables HI and ITT1 Both the widths and the over- 
all yields of pions increase as a function of beam energy. 
However, there is a more significant increase in the ob- 
served pion yields between 2 and 4 AGeV than between 
4 and 6 AGeV or 6 and 8 AGeV. 



Efceom Q^- (7T~) W^- 

2 AGeV 21.3 ± 0.1 ± 1.3 36.1 ± 0.3 ± 2.0 0.675 ± 0.006 ± 0.005 
4 AGeV 39.0 ± 0.1 ± 2.1 76.0 ± 0.2 ± 4.2 0.780 ± 0.003 ± 0.002 
6 AGeV 50.8 ± 0.1 ± 2.7 104.0 ± 0.2 ± 5.4 0.817 ± 0.002 ± 0.001 
8 AGeV 61.1 ± 0.1 ± 3.3 130.7 ± 0.4 ± 7.9 0.854 ± 0.003 ± 0.008 

TABLE I: dN/dy fit parameters for negative pions at each 
beam energy with statistical and systematic uncertainties re- 
ported separately. Q„.- is the amplitude of the distribution 
a t y-ycM = 0, (tt~) is the 4-n yield (total area) and w w - is 
the width. 



IV. DISCUSSION 

In this section, some of the characteristics of the ob- 
served pion spectra obtained from E895 are discussed. 
The pion spectral shapes are compared to the predictions 
of a microscopic transport model, the dN/dy distribu- 
tions are evaluated in the context of collective dynamics, 
and the overall yields are used to infer the initial state 
entropy density obtained in these collisions. 

It is interesting to note the asymmetry between the 
positive and negative pion yields at each beam energy. 
Although the ratio of negative to positive pions decreases 
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FIG. 6: Pion transverse mass spectra at mid-rapidity for 2, 
4, 6, and 8 AGeV (uppermost curves on Figs. 0]and|^. Two- 
slope thermal fits are shown superposed on the data along 
with the contributions to the total from the two separate in- 
verse slope terms in Eq. 10 ■ 



FIG. 8: Rapidity dependence of the fit parameters from 4 
AGeV n~ mt — mo spectra fits with a two-slope model. See 
caption of Fig. [7| for details. 
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FIG. 7: Rapidity dependence of the fit parameters from 2 
AGeV tv~ mt — mo spectra fits with a two-slope model. The 
leftmost panels show the inverse slope parameters. Ti (indi- 
cated by closed triangles) dominates the low-(m t — mo) range 
of the spectra, whereas T2, (indicated by closed squares) dom- 
inates the high-(m t — mo) end of the spectra. The middle 
panels show the two dN/dy fit parameters, (dN/dyi in open 
triangles and dN/dy2 in open squares) and the rightmost pan- 
els show the total dN/dy, which is the sum of dN/dy 1 and 
dN/dy2. The top row reports the fit results when all four 
parameters are allowed to be independent in the fit. The bot- 
tom row shows the results for the remaining three parameters 
when the high-(mt — mo) inverse slope is fixed to the Gaussian 
estimation. Note that the total dN/dy is fairly insensitive to 
the interplay among the parameters. 
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FIG. 9: Rapidity dependence of the fit parameters from 6 
AGeV 7r~ mt — mo spectra fits with a two-slope model. See 
caption of Fig.0for details. 



2 AGeV 11.5 ± 0.3 ± 1.2 19.2 ± 1.3 



+2.0 
-1.0 



0.668 ± 0.053 ± 0.011 



4 AGeV 27.7 ± 0.3 46.3 ± 0.8 ttl 0.667 ± 0.015 ± 0.025 
6 AGeV 38.4 ± 0.3 ±jj 75.7 ± 1.1 t 3 2 .t °- 787 ± °- 014 ± °- 016 
8 AGeV 46.2 ± 0.4 ±l\ 95.9 ± 1.1 tH 0.828 ± 0.013 ± 0.005 

TABLE II: dN/dy fit parameters for positive pions at each 
beam energy with statistical and systematic uncertainties re- 
ported separately. is the amplitude of the distribution 
at y-ycAf = 0, (tv + ) is the 4-n yield (total area) and w„.+ is 
the width. 
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FIG. 10: Rapidity dependence of the fit parameters from 8 
AGeV 7T _ mt — mo spectra fits with a two-slope model. See 
caption of Fig. [7| for details. 
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FIG. 11: Rapidity dependence of the fit parameters for each 
beam energy from the 7r + mt — mo spectra fits with a two- 
slope model. (See caption of Fig. 0for more details.) The 
high-(mt — mo) inverse slope parameters have been fixed to 
the 7T~ parameterized values. 



over the studied beam energy range, from 1.88 at 2 AGeV 
to 1.36 at 8 AGeV (see Tables III and Hill . it does not reach 
the asymptotic value of ~ 1.0 observed at the top CERN 
SPS energy [3 S3- The neutron excess in Au+Au col- 
lisions (118+118 neutrons compared to 79+79 protons) 
combined with the pion branching ratios suggests that 
1.95 negative pions for every positive pion will be pro- 



FIG. 12: Yield per event of 7r _ (left-hand panel) and 7r + 
(right-hand panel) in 0-5% central Au+Au collisions at 2, 4, 
6, and 8 AGeV integrated over m t — mo as a function of rapid- 
ity. A Gaussian parameterization is applied to extract the An 
yields. Reported uncertainties are statistical plus systematic. 



duced . If the isospin fractions are folded with the ob- 
served cross-sections for NN — > NN-7T from experimental 
measurements |40j . the expected ratio of (7r~):(7r + ) at 1 
AGeV is ~ 1.91:1. The 2 AGeV data, with ^/IWn not far 
from the A(1232) production threshold, are quite near, 
though a little lower than this predicted ratio. As the 
energy increases, the number of directly produced pion 
pairs (7r - 7r + ) is expected to lower the ratio, asymptoti- 
cally approaching 1.0, which is the trend we observe. At 
8 AGeV the negative pion excess is only approximately 
36%, compared to 88% at 2 AGeV. 



A. RQMD Comparisons 

RQMD (Relativistic Quantum Molecular Dynamics) 
v2.3 |4l| is a microscopic transport model which attempts 
to simulate heavy ion collisions by propagating individ- 
ually all particles through the six dimensions of phase 
space in the fireball. Interaction probabilities are ap- 
proximated by using published interaction cross-sections 
of free hadrons and the relative phase space proximity 
of pairs of particles at each time step of the reaction. 
Inelastic collisions may produce new particles, such as 
pions, which are also propagated through phase space 
along with the nucleons. The reaction ends when the 
phase space density reaches a low enough threshold such 
that the probability of further interactions is small - the 
freeze-out point. This model has been reasonably suc- 
cessful in describing many final state observables experi- 
mentally measured in the beam energy range studied for 
this analysis [42| . 

In this and other cascade models, final-state particle 
distributions are frozen at the end of the reaction. Post 
freeze-out effects, such as the Coulomb interaction of the 
pions with the nuclear fireball, which will be discussed 
in detail using results from this analysis elsewhere |38| . 
are not included. However, RQMD combined with an 



afterburner to simulate final state Coulomb interactions 
|43| and to permit the decay of residual resonances has 
been successful at describing asymmetries in observed 
pion yields at low p t . 

RQMD version 2.3, with the nucleon mean field set- 
ting turned on, was used to generate pion distributions 
from central (b < 3 fm) Au+Au collisions for compari- 
son with the data obtained by E895 for this paper. One 
of the output parameters from RQMD records the na- 
ture of the last collision, "lastel" , of each particle before 
freeze-out. Particles whose lastel was a thermal rescat- 
tering are labelled "thermal" . Particles whose lastel was 
the decay of a resonance have the parent particle listed 
explicitly. However, high mass resonances which are not 
Deltas or vector mesons are combined in a single cate- 
gory labelled "himass" in the user's notes for the code. 
Particles whose lastel was String or Rope fragmentation 
are indicated separately. Figs. 1131 and ITU show the na- 
ture of the last collision for ir~ and tt + for a centrality 
selection of b < 3 fm at 2 AGeV and 8 AGeV, respec- 
tively. The distributions are normalized to the per event 
yield of pions. In each case, the Delta resonances are the 
largest single contributor to the last pion interactions be- 
fore freeze-out. 

The dependence of the pion spectral shapes on feed- 
down from Delta resonances is evident in Fig. 1151 which 
shows the m t — too spectra of pions near mid-rapidity 
(|y| < 0.3) from RQMD. Pions whose last collision was 
a Delta are plotted separately from the pions with all 
other last interactions. The spectral shapes are strongly 
affected by a Delta lastel. Plotted in Fig. are the 
ratios of the pion yields from Deltas and all other last 
interactions to the total yields as a function of mt — too. 
In all cases the Delta pions contribute more significantly 
to the total yield at low m f — too than at high m t — too, 
but the influence of the Delta contribution to the spectra 
diminishes as the beam energy increases. The total mid- 
rapidity pion spectra are fit with the same model as was 
used to describe the data, Eq. ©, with the fit parame- 
ters shown in Tables IIIII and IIVI At both high and low 
TO* — mo, the RQMD ir~ and ir + appear to have common 
inverse slope parameters and there is no obvious evolu- 
tion with beam energy. The average value at low m t — too 
(Ti) is —72 MeV, while the high m t — to (T 2 ) average 
is -140 MeV. The magnitude of the low m t - m RQMD 
inverse slope parameters is much larger than what is ob- 
served in the data at all beam energies. At high m t — too, 
the RQMD inverse slope parameters are much closer to 
the observed values, except at 2 AGeV, where the data 
show T 2 ~ 120 MeV at mid-rapidity. 

The rapidity density distributions of charged pions pre- 
dicted by RQMD are shown in Fig.ll7lalongside the E895 
data from Figs. 1711 II The contributions to the total 
dN/dy from Delta resonance decay pions and all other 
pions are also shown. As the beam energy increases, the 
contribution to the total from Delta resonance decay pi- 
ons becomes less important, which is consistent with the 
trends seen in the data. 
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TABLE IV: Mid-rapidity RQMD tt + fit parameters, Eq. 0. 



10.2 


73 


1.6 


134 


14.4 


72 


8.5 


131 


15.5 


64 


17.7 


131 


21.4 


70 


21.3 


144 



o 
s 
a 

T3 

2: 
-a 



100 
10 

1 

io-t 

10~ 3 
10 

1 

10t 

lo- 
rn 1 

"3 



2 AGeV 71" 



2 AGeV 7C + 



.a « 

b B 

& £ 

H - 



RQMD v2.3: Nature of last collision before freezeout 



FIG. 13: RQMD nature of pion last collision, "lastel", before 
freeze-out at 2 AGeV for impact parameter, b < 3 fm. 
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FIG. 14: RQMD nature of pion last collision, "lastel", before 
freeze-out at 8 AGeV for impact parameter, b < 3 fm. 
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FIG. 15: RQMD charged pion spectra at mid-rapidity for 
collisions with impact parameter, b < 3 fm. The contribution 
to the spectra from pions whose last interaction was a Delta 
decay and the sum of all other processes are shown separately. 
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FIG. 16: Ratios of pion yield for different production mech- 
anisms to the total yield as a function of mt — mo at mid- 
rapidity from RQMD collisions with impact parameter, b < 
3 fm. The contribution to the spectra from pions whose last 
interaction was a Delta decay and from the sum of all other 
processes are shown separately. 



B. Longitudinal flow 

An important question in heavy ion collisions is the de- 
gree of collectivity of the produced particles, which can 
arise from the build-up of pressure in the hot, dense col- 
lision zone. RQMD has no such hydrodynamic effect ex- 
plicitly included. However, individual particle thermal 
rescattering, the microscopic analog to pressure, can also 
drive collective flow, and has been observed in RQMD 
calculations at RHIC energies 0]. Rapidity broaden- 
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FIG. 17: Rapidity density of charged pions from RQMD cal- 
culations for collisions with impact parameter, b < 3 fm. The 
contribution from Delta resonance decay pions and the sum of 
all other processes are also plotted separately to show the rel- 
ative contribution to the pion yields from Delta decays. Also 
plotted are the data from Figs. 171111 



ing along the beam axis has been observed in heavy ion 
collisions at many beam en erg ies for a wide range of sys- 
tems and particle species [22, |23|, 0, EH • This broad- 
ening has been interpreted as arising from the collec- 
tive motion of the system after collision in the longitu- 
dinal direction (longitudinal flow). The observed rapid- 
ity distributions are compared with the expectation for 
a stationary thermal source, and with a longitudinally 
boost-invariant superposition of multiple boosted indi- 
vidual isotropic, locally thermalized sources in a given 
rapidity interval. Each locally thermalized source is mod- 
elled by the m t -integrated Maxwell-Boltzmann distribu- 
tion, Eq. i with the rapidity dependence of the energy, 
E = m t cosh(y) explicitly included, and T is true tem- 
perature: 

dN th 



\ e (- ! T coshy) 

(5) 

The distributions are integrated over source element 
rapidity to extract the maximum longitudinal flow, rj max 

dy 

0L 



, - „ . . .0 , ill lib Zj £j 

-r- 60 = BT 3 (— + + —Tr- 

ay T z T cosh y cosh y 



, dN th , 
dr)——[y-ri) 
dy 



(6) 



tanh(?7. 



0- 



where T) max — - rj m i n , from symmetry about the center 
of mass, and (3l is the maximum longitudinal velocity in 
units of c. An average longitudinal flow velocity can be 
defined as (/?£,) = tanh(r) max /2). 

Ref. HH presents longitudinal flow velocities extracted 
from the proton rapidity densities measured by E895 for 
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the same event selection as the pions presented in this 
analysis and compares the results to values extracted 
from a wide array of experiments over the beam energy 
range from 1 to 160 AGeV. Since the protons are present 
before the collision, at higher beam energies they may 
be strongly affected by nuclear transparency, which can 
also broaden the rapidity distributions. Therefore, the 
protons alone cannot be used to determine the absolute 
magnitude of the collective motion, if it is present. 

In order to determine the degree of collectivity, it is 
important to compare multiple particle species from the 
same collisions simultaneously. This has been done with 
central Si+Al collisions at 14.6 AGeV for protons, pi- 
ons, kaons and lambda hyperons (Z^ and in central S+S 
collisions at 200 AGeV for pions, kaons and lambda hy- 
perons [23|. A common collective flow velocity was able 
to reasonably describe all of the observed rapidity dis- 
tributions, except the protons in S+S collisions at 200 
AGeV, where apparent nuclear transparency is more pro- 
nounced. A similar simultaneous description is possible 
here, by combining the proton rapidity distributions from 
|23|| with the present pion rapidity distributions. 

The rapidity densities of pions and the protons from 
Ref. 2.3] are shown in Figure El Stationary thermal 
source emission functions, the sum of two of Eq. (J5J for 
the pions using the two inverse slope parameters from 
the pion transverse mass spectra fits at mid-rapidity, are 
shown as dashed lines. We have not measured the true 
system temperature, T, since Ti and T2 parametrize 
the temperature and known radial flow, resonance and 
Coulomb effects in the pion spectra. Consequently, T2 
overestimates the system temperature and Ti may un- 
derestimate (due to competing effects of radial flow, res- 
onance feed-down and the Coulomb interaction with the 
fireball). Therefore the plotted distributions (dashed 
lines) are probably wider than the true thermal distribu- 
tions would be, and yet are still too narrow to reproduce 
the observed rapidity spectra. The inclusion of longitudi- 
nal flow remedies this. For the protons, a single Eq. (J5J 
is used for the thermal rapidity distribution, as in Ref. 

M 

The solid curves are the emission functions including 
longitudinal flow, Eq. JSJ, with the velocities fixed to 
the values extracted from the protons ((/3l) = 0.28, 0.42, 
0.48, 0.50 at 2, 4, 6 and 8 AGeV, respectively). The longi- 
tudinally expanding source clearly better reproduces the 
measured pion distributions at all beam energies than do 
the dashed curves in Fig. 1181 Consistency among parti- 
cles of different masses supports a hydrodynamical inter- 
pretation of the rapidity density broadening; the system 
is expanding like a fluid with a common longitudinal flow 
velocity. 



C. Entropy production 

It was sug gest ed many years ago by Fermi |46j and 
later Landau |47| that pion production may be used to 
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FIG. 18: 7T~, 7r + and proton rapidity density distributions 
at 2,4,6, and 8 AGeV. The expectation for isotropic emission 
from a stationary thermal source is indicated with dashed 
lines, whereas the solid curves represent the form including 
longitudinal flow from the proton rj m ax values given in Ref. 
I2J In all cases the thermal model is too narrow. The proton 
data and the longitudinal expansion velocity come from Ref. 
0. 



estimate the amount of entropy produced in high energy 
particle collisions. Later, Van Hove ji^] extended this 
idea to heavy ion collisions and proposed that this may 
be a way to distinguish events in which a Quark Gluon 
Plasma (QGP) is formed. In a QGP the color degrees of 
freedom of the liberated partons introduce a significant 
number of new energy states unavailable in a hadron gas. 
By studying pion production over a broad range of colli- 
sion systems and energies, discontinuities in the observed 
multiplicities might indicate the onset of QGP formation. 

In 1995, Gazdzicki ^| took the available data from 
heavy ion collisions and showed that there is an increase 
in the observed entropy produced at the SPS (NA35 Ex- 
periment with S+S collisions at 200 AGeV 50] ) com- 
pared with AGS energies. The low energy heavy ion data 
follow the trend for p+p collisions, while at the SPS there 
is an apparent factor of 3 increase in the effective num- 
ber of degrees of freedom |5lj| . The entropy production 
analysis of Pb+Pb collisions at 158 AGeV (NA49) H2 
supported this observation. 

This model assumes that the entropy is produced at 
the early stage of the collision when the incident matter 
is in a highly excited state. The thermalized, strongly 
interacting matter is assumed to expand adiabatically to 
the freeze-out point, preserving the early stage entropy. 

Since the majority of produced particles are pions, to 
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first order, the mean pion (boson) multiplicity should 
be nearly proportional to the entropy. The ratio of the 
mean pion multiplicity to the mean number of participat- 
ing nucleons, (it) / (N part } , provides a simple estimate of 
the entropy density. (N part ) for a nucleus-nucleus colli- 
sion can be estimated using a Glauber model calculation 
of the mean free path of the nucleons through the nu- 
clei as they collide at a given impact parameter. For 
the present analysis, (Np) was estimated using RQMD. 
The nucleons are distributed according to a Woods-Saxon 
nuclear density profile and the impact parameter of the 
simulated collision is used with the nucleon-nucleon in- 
teraction cross-sections [53| to determine the number of 
participants. The top 5% of collisions (determined by 
integrating the (Np) distribution from a set of minimum 
bias RQMD events at each beam energy) correspond to 
(N P ) = 364, 366, 365, 363 at 2,4,6, and 8 AGeV, respec- 
tively. The estimated uncertainty on these values is ± 5 
participants. 

Following Ref . [24[ , the entropy densities for each beam 
energy, here approximated as (it) / (N part ) , with (n) = 
1.5((tt~) + (tt + )) to account for the neutral pions, are 
plotted in Fig. ^| as a function of the Fermi energy vari- 
able, F = (vg^Zzggw) 3 / l_ E895 data are indicated by 

stars and the values are tabulated in Table E| NA49 
results for 40, 80 and 158 AGeV Pb+Pb collisions from 
the CERN SPS were obtained from Ref. . The number 
of participants from 0, calculated using the FritiofjH^I 
model, are also listed in Table IVl 

The linear dependence of the entropy per participant 
nucleon as a function of F in the proton-proton(anti- 
proton) data is not evident for the full range of the heavy 
ion collision data. At and below AGS energies, the heavy 
ion data lie below the p+p data, and appear to be approx- 
imately linear with F. In Ref. [24j, the SPS results com- 
bined with RHIC results at much higher energies from the 
PHOBOS Collaboration show a linear trend with a slope 
that is approximately 1.3 times larger than at the lower 
energies. There appears to be a transition in the region 
between the AGS and top SPS energies. The third-order 
polynomial fit to the heavy ion data shown on Fig. 1191 
may indicate that a smooth trend with increasing F can 
accurately describe the excitation function without the 
need for a discontinuous jump, such as one might expect 
from a first-order phase transition. Two more runs at 
the SPS with beam energies of 20 AGeV and 30 AGeV 
may be able to improve the resolution in this important 
transition region. 



V. SUMMARY 

Transverse mass and rapidity spectra of charged pi- 
ons in 2-8 AGeV 0-5% central Au+Au collisions have 
been measured by the E895 experiment. The transverse 
mass spectra exhibit a low-pt enhancement which can be 
largely ascribed to the feed-down from late stage reso- 



Vsaw (GeV) F (GeV 1/2 ) 


{Np) 


(n)/{N p ) 


2.630 


0.644 


364 


0.2279 ± 0.0159 tomei 


3.279 


0.965 


366 


0.5012 ± 0.0111 


3.838 


1.190 


365 


0.7385 ± 0.0148 


4.289 


1.351 


363 


0.9364 ± 0.0170 ±0.0579 


8.830 


2.452 


349 


2.6433 ± 0.0858 


12.280 


3.099 


349 


3.9542 ± 0.0869 


17.260 


3.821 


362 


5.2127 ± 0.1823 



TABLE V: Tabulated mean number of pions per participant 
and Fermi energy variable for 2, 4, 6, and 8 AGeV Au+Au 
collisions and 40, 80, and 158 AGeV Pb+Pb collisions from 
Ref. 24]. Statistical and systematic uncertainties arc reported 
separately. 
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FIG. 19: Mean pion multiplicity per participant vs. F, the 
Fermi energy variable. The solid symbols represent data from 
heavy ion collisions, whereas the open symbols come from 
p+p interactions. The E895 data points are indicated by 
stars. The third-order polynomial fit to the heavy ion data 
show that a smooth trend with F is possible, although the 
low energy and high energy data (including RHIC results) can 
be well-described by separate linear parameterizations whose 
slopes differ by a factor of approximately 1.3 24]. 



nance decays. Differences in the tt + and 7r~ spectra at 
low rrit — mo are not reproduced by RQMD, which does 
not include final state interactions such as the Coulomb 
interaction with the nuclear source. The inverse slope 
parameters increase as a function of beam energy and 
appear to be charge independent at high m t — too. The 
measured rapidity distributions show excellent forward- 
backward rapidity symmetry and are well described by a 
model which includes collective longitudinal flow, with a 
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velocity that is common to both pions and protons emit- 
ted in these collisions. The 4-7T yields of pions, obtained by 
integrating the rapidity distributions, have been used to 
infer an initial state entropy which increases with beam 
energy and is consistent with a smooth non-linear trend 
as a function of F from the 2 AGeV Au+Au collisions to 
158 AGeV Pb+Pb collisions at the SPS. 
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